Dynamic phase transition in a kinetically constrained model for traffic 
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We address the emergence of traffic jams using dynamic correlation functions in analogy to kinet- 
ically constrained models for glasses. In kinetically constrained models, the formation of glass be- 
comes a true (singular) phase transition in the limit T — > 0. Similarly, using the Nagel-Schreckenberg 
model to simulate traffic flow, we find that the emergence of jammed traffic acquires the signature 
of a sharp transition in the deterministic limit p —¥ 1, corresponding to overcautious driving. We 
identify for the first time a true dynamical critical point marking the onset of coexistence between 
free flowing and jammed traffic. We find diverging correlations analogous to those at a critical point 
of thermodynamic phase transitions. 
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Dynamic arrest, the sudden slow-down of dynamical 
systems with increasing density or interaction potential, 
is a central phenomenon in complex systems across bi- 
ology, geology, material science, transport and traffic. 
While the dynamic arrest is important for material sta- 
bility and memory, it is rather detrimental in traffic or 
transport, where congestion freezes any motion. A well- 
known example of dynamical arrest is the sharp increase 
of the viscosity of glass forming liquids at the glass tran- 
sition. Many more systems exhibit dynamic arrest phe- 
nomena that are seemingly similar. Examples include 
dynamics in the crowded environment of cells or human 
dynamics in escape situations or traffic flows. 

While the nature of this transition is not clear, a pos- 
sible idea is that of a dynamic phase transition, in which 
the dynamic slow down is considered to be analogous to 
equilibrium phase transitions with their singularities in 
thermodynamic quantities. An important question then 
concerns the universality of the dynamic arrest. While 
dynamic arrest phenomena across different systems seem 
to be related, it is not clear to which extent they really 
are. 

Important insight into the dynamic arrest of glasses 
comes from kinetically constrained models (KCMs), a 
class of discrete models with stochastic dynamics that 
are used to describe the glassy behavior and increas- 
ing relaxation time scales in supercooled liquids As 
the defining ingredient, KCMs have a kinetic constraint 
that allows local activity only if a local condition is met. 
These models provide some evidence that indeed the dy- 
namic slow-down is the manifestation of a dynamic crit- 
ical point @, 0] in the limit T -> 0. 

A constraint analogous to those in KCMs exists in traf- 
fic flow: cars can accelerate only if the distance to the car 
in front is sufficiently large. A well-studied model that 
incorporates a number of basic dynamical properties of 
real traffic is the Nagel-Schreckenberg (NS) model a 
lattice-gas-like model with discrete position, time, and 
velocity. It allows for stochastic fluctuations in the ve- 
locities of the individual cars, controlled by the proba- 
bility p that reflects the drivers' individual freedom to 



adjust their speed. The NS model describes the forma- 
tion of traffic jams, and there has been much discussion 
regarding the (non)cxistcncc of a sharp phase transition 
between free- flowing traffic and traffic with jams 

We reserve the word dynamical phase transitions here 
to denote qualitative changes of behavior, accompanied 
with singularities in various expectation values in the 
stationary state. This is analogous to the singularities 
of observables at thermodynamic phase transitions, with 
the equilibrium state replaced by the stationary state. 
While different types of traffic flow are often referred to 
as "phases" Ref. p, 0], the traffic jams have a typical 
size, and there is no divergence of correlation lengths [|| , 
as one expects to be at a phase transition. Consequently, 
the evidence of dynamical criticality and the connection 
to other dynamically arrested systems remains elusive. 
However, a sharp transition without diverging correla- 
tions was found in the limit p — ¥ jjj ■ 

In this Letter, we show that in the deterministic limit 
p — > 1, the simple one-dimensional NS model 0] of traf- 
fic flow does exhibit a true dynamical phase transition. 
We apply concepts normally used for glasses, such as 
dynamic correlation functions and susceptibilities Q to 
show that in the deterministic limit, the maximum corre- 
lation length and time diverge, analogous to the dynamic 
transition in KCMs at T — > 0. These relatively sim- 
ple stochastic dynamical systems allow for a direct com- 
parison of dynamic arrest in different dimensions. We 
identify observables in traffic which behave similarly and 
which relate to each other in the same way as in glasses. 
This signals some degree of universality in dynamic ar- 
rest. 

The NS model simulates traffic flow on discrete space 
and time. A fixed number of cars with average density p 
per lattice site are positioned on a one-dimensional lattice 
with periodic boundary conditions. The cars have inte- 
ger velocities vi between and some maximum velocity 
» max . The dynamics is given by the following update 
rules applied in parallel to all cars: All cars i with veloc- 
ity Vi < Umax accelerate by 1. Next, the kinetic constraint 
is applied, so that, if the distance of car i to the next car 
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FIG. 1: Phase diagram of the Nagel-Schreckenberg model 
showing (v) for the case of t) max = 2 (see gray scale on the 
right). The dashed line marks the transition between freely 
flowing traffic and traffic with jams found in Eq. |T}. 



di < Vi, then car i decelerates to di- With probability 
p, a car reduces its velocity by 1. This parameter is the 
only source of stochasticity in the system and controls 
the fluctuations of the velocities of cars. Finally, all cars 
move along the road by vi lattice sites. 

We give an overview of the average velocity (v) as a 
function of density p and probability p in Fig. [TJ For sim- 
plicity, we have taken u max = 2. At low density, traffic 
can flow freely, and (v) = v f = i' max — p. With increas- 
ing density, cars interact and decelerate according to the 
kinetic constraint, leading to the formation of jams and 
a decrease in (v). We define the transition density ptra 
from the balance of the outflow and inflow rates of a jam, 
as required for its stability. The outflow rate of a contin- 
uous sequence of jammed cars is 1 — p, corresponding to 
the probability of acceleration of the car at the head of 
the jam. Because cars approach the rear of a jam with av- 
erage velocity Vf = v max ~P, and the rear of the jam itself 
travels backwards at a speed equal to the inflow/outflow 
rate, this yields the transition density 
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Consequently, p tla , oc 1 — p for small 1 — p. 

For glasses and KCMs 0, 11, 12 1, it has been shown 
that the transition becomes singular at T = 0, where the 
dynamics become deterministic. The question is then 
whether a similar singular transition exists in a simple 
one-dimensional model of traffic flow. To explore this 
analogy, we relate the stochasticity parameter p in the 
NS model to the temperature T of spin glasses. The case 
T = 0, where the dynamics of glasses freezes entirely, 
corresponds to the case p = 1, where cars always decel- 
erate, and traffic flow arrests. In the limit of T — > 
and p — > 1, the systems become deterministic. For real 



traffic flow, the case p — >• 1 is rather exceptional, but 
can be approached in the case of bad weather conditions, 
where drivers act carefully, and are likely to overreact 
with braking. 

We investigate the onset of jamming at p ~ p tra , by 
applying dynamic correlation functions that have been 
much used to investigate the dynamic arrest of glasses. 
A characteristic property of glasses is their dynamic het- 
erogeneity. Dynamically active regions separate from dy- 
namically less active regions in space and time, leading 
to increasing dynamic heterogeneity of the system. This 
dynamic heterogeneity can be quantified by the dynamic 
susceptibility 0, El, EH ■ We define the dynamic correla- 
tion function of traffic flow using 



G 4 (i,t) = (c(i;t)c(0;t)) - (c(0;i)> : 



(2) 



where we take the mobility c(i;t) of car i as c(i;t) — 
(l/{t + 1)) X)t'=o its average velocity during the 

time interval [0, t] . We find that near p = 1, dynamic cor- 
relations become indeed increasingly long-ranged when 
the density approaches ptr&, the onset of the jammed 
regime. To investigate whether there is a true diverg- 
ing length scale, we define the dynamic susceptibility 
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that measures the number of cars that move coopera- 
tively on the time scale t. The dynamic susceptibility %4 
indicates the size of regions of correlated mobility, and 
has been much used to measure dynamic heterogeneity 
in glasses and granular materials [15l-ll8|. 

In glasses, maximum cooperative motion arises at in- 
termediate time scales, at which the particles escape their 
cages. In traffic, cars cannot escape the crowding of 
their environment independently from each other, and 
the maximum dynamic susceptibility arises at the short- 
est time interval, t = 1, sec inset of Fig. [2^. We use the 
maximum value of %4 to explore dynamic correlations for 
a range of p- values in Fig. [2ta). Indeed increasing max- 
ima develop at p ~ pt ra , as p approaches unity, indicating 
increasing dynamic correlations. The divergence of the 
dynamic susceptibility is clearly seen in Fig. [2Jd, where 
we plot the maximum value of x<t as a function of 1 — p. 
The figure shows data for various v max ; in all cases, the 
maximum Xi grows as a power law X4,max = (1 — P)~ u , in- 
dicating that the number of cars that move cooperatively 
diverges. The divergence of the dynamic susceptibility 
indicates that traffic flow in the NS model becomes truly 
critical in the limit p — > 1. In this model, the exponent 
of the power law appears to increase weakly with u max 
changing from v = 0.53 to v > 0.70. We also note that 
no maximum of %4 is observed for p — > 0, indicating the 
absence of dynamic correlations, and the very different 
nature of this limit. The criticality for p — )■ 1, however, 
is analogous to the one observed in KCMs at T — > 0, and 
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FIG. 2: Dynamic susceptibility of traffic flow in the NS 
model, (a) The value of X4(0) as a function of rescaled den- 
sity for a range of probabilities p at v msx = 2. The peak 
sharpens markedly as p — > 1. Inset: X*(f) as a function of 
time shows that the largest value occurs at t — 0. (b) Maxi- 
mum value of \4,(t = 0) as a function of the rescaled density, 
plotted vs. I — p, for a set of v max (symbols explained in (c)). 
Power-law behavior (dashed lines) indicates the divergence of 
the dynamic susceptibility on approach of the critical point 
p — 1. (c) Density of maximum dynamic susceptibility as a 
function of 1 — p for various values of D max . The position of 
the maximum is well described by the limiting behavior of 
(HJ), Ptra = (1 — p)/(wmax — 1) indicated by the dashed lines. 
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FIG. 3: Residence time of cars in jams as a function of 
reduced density for various values of p at i> max = 2. As p — > 1, 
the correlation time becomes more sharply peaked. Inset: 
Maximum residence time as a function of 1 — p shows power- 
law divergence of the correlation time for p — > 1 . 



indicates that the deterministic limits p — > 1 and T — > 0, 
are dynamical critical points of the systems. The diver- 
gence occurs at the onset of the jammed regime; to show 
this, in Fig. [5{c), we compare the location of the maxi- 
mum of %4 (symbols) with the limiting (p — > 1) behavior 
of ptra according to Eq. [T] (dashes lines). 

Further evidence of critical behavior comes from mea- 
surement of the correlation time scale. To estimate the 
typical time scale, we make use of a quantity similar to 
%4, where we interchange time and car index in the defi- 
nition of c(i; t) and in the sum appearing in ([3]), to obtain 
$4(1). The temporal susceptibility "$4 indicates the corre- 
lation time scale of the system, and measures the typical 
residence time of a car in a jam. Wc plot this corre- 
lation time as a function of reduced density in Fig. [3] 
A strong increase of the maximum of $4 suggests that 
in addition to the divergence of the correlation length, 
there is also a divergence of the correlation time scale. 
This is confirmed by plotting the maximum values of $4 
as a function of (1 — p) in the inset. Similar to the spa- 
tial correlations, the correlation time scale diverges as 
a power law ^4 ;ma x = (1 — as p — > 1, confirming 
that the system behaves critically along the time dimen- 
sion. We determine the exponent to be /1 ~ 1.5. For real 
traffic, such diverging correlation time can have fatal con- 
sequences, as it indicates diverging persistence times of 
traffic jams. There are historical examples of traffic jams 
lasting several days (see, for instance, Ref. [lj|), under 
particular weather conditions. 

We thus find a dynamical critical point with diverg- 
ing length and time scales. We note that the condition 
of densities large enough to cause jams, but not large 
enough to jam the whole road over its entire length, is 
analogous to the coexistence of two phases. Indeed, in 
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FIG. 4: Order parameter M as a function of reduced density 
for various values of p at u max = 2. As p — !> 1, the transition 
at p — 1 — p becomes singular. The dashed line shows the 
theoretical prediction M = (p — ptm) / P- 



the limit p —> 1, jams coalesce in time to form a single 
jammed phase coexisting with a single free flowing phase. 
Taking into account the divergence of the quantites \i 
and z?4 on entry of the coexistence regime, the situation 
appears similar to a thermodynamical critical point ter- 
minating a two-phase regime. 

To explore this analogy further, we define a dynamical 
order parameter 



M = 



Vf - (v) 



(D 



in terms of the average velocity (v) and the mean veloc- 
ity in free flow Vf = l> max — P- We show M as a function 
of the rescaled density in Fig. It is continuous at the 
transition p = ptra, but exhibits an increasingly sharp 
kink as p — » 1, indicating a singular point in the limit 
p — > 1. If wc assume simple coexistence in the two-phase 
regime, we can predict the function M(p) = [p — ptm)/p- 
Indeed as p — >• 1, there is strong evidence that the data 
converge to this simple function, supporting our picture 
of jam and free flowing traffic as coexisting phases. The 
functional dependence of this order parameter has an ex- 
ponent /8 = 1, corresponding to a Bose condensate, and 
to condensates found in typical zero range models [20|. 

We have shown that a simple ID model for traffic flow 
exhibits signature of dynamic arrest analogous to that 
of kinetically constrained models for glasses. Exploit- 
ing this analogy, we have identified and analyzed a non- 
equilibrium phase transition in the Nagel-Schrcckcnberg 
model in the limit p — > 1 and p — > 0, between a free 
flowing state and coexistence of jammed and free flowing 
traffic. The hallmark of this transition is the divergence 
of both correlation length and time scales, giving it a 
second order character. Defining a proper dynamic order 



parameter, we have shown that this transition is analo- 
gous to a thermodynamic phase transition in which the 
coexistence regime is entered through the critical point. 
Finally, the direct analogy to KCMs of glasses points out 
the universality of dynamic arrest phenomena in systems 
of different dimensionality. 
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